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ABSTRACT 

We present for the first time orbital phase-resolved spectra of an intermediate po- 
lar (IP), EX Hya, together with the spin phase- resolved spectra during two different 
epochs using the X-ray Multi-Mirror Mission (XMM-Newton ), European Photon 
Imaging Camera (pn instrument). We find that the source at the two epochs has the 
same X-ray luminosity of ~ 6.5 x 10 31 erg s _1 . We detect spectral variations between 
the 2000 and 2003 observations of the source. We fitted the spectrum using a neutral 
hydrogen absorption model with or without covering fraction together with Gaussians 
for emission lines, two collisional equilibrium plasma emission models (MEKAL) and a 
cooling-flow plasma emission model (VMCFLOW). We find that two of the three emis- 
sion components (>T=0.6-0.8 keV and fcT=1.3-I.7 kcV) fitted by the MEKAL models 
are almost constant over the spin and orbital phases and also over the two different 
epochs with the normalisation varying directly proportional to the flux when the data 
are folded according to the orbital and spin phase indicating that the slight variation 
may be due to occultation. The emission modeled by the VMCFLOW changes over 
the spin and orbital phases and the 2000 and 2003 observations reveal two different 
ranges of temperatures (3-33 and 8-61 keV respectively) that model the shock zone in 
the accretion column/s. The ratios of the spin maximum to minimum and the orbital 
maximum to minimum spectra along with the increase in the plasma temperatures 
indicate that the spectrum gets harder in the minimum phases of both orbital and 
spin periods. In the 2003 observation, a 6.4 keV fluorescent Fe emission line is present 
at the orbital minima in a range of phases from 0.9 to 1.3 and it is absent otherwise. 
This indicates that there is reflection from the disc most likely from a large bulge at 
the accretion impact zone. 

Key words: binaries: close - Stars: individual: EX Hya - Intermediate Polars, cata- 
clysmic variables - stars:rotation - white dwarfs - X-rays: stars 



1 INTRODUCTION 

Cataclysmic variables (CVs) are compact binaries hosting a 
white dwarf (WD) primary star accreting material from the 
Hydrogen rich companion. The accretion occurs through an 
accretion disc in cases where the magnetic field of the WD 
is weak ( B < 0.01 MG) and such systems are referred to 
as nonmagnetic CVs. The magnetic CVs (mCVs) consti- 
tute about 25% of the entire population of CVs. There are 
two subclasses depending on magnetic field strength. Polars 
(AM Her type systems) have B > 10 MG where this strong 
field causes the accretion flow to directly channel onto the 
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magnetic pole of the WD inhibiting the creation of an ac- 
cretion disk (see Cropper 1990; Warner 2003). Polars show 
strong circular and linear polarization modulated at P or b 
and the WD rotation P sp i n is synchronized with the binary 
orbit. Only about 1% of the Polars are found to show asyn- 
chronicity. The second class of mCVs are Intermediate Po- 
lars (DQ Her type systems) that have a weaker field strength 
compared with Polars of about 1-10 MG. In this case, the 
accretion takes place via a truncated disk to the magnetic 
poles through accretion curtains. Near the surface of the 
WD, the accreting material forms a strong shock where the 
post shock region heats up to 10-20 keV and then cools 
via thermal Bremsstrahlung (Patterson 1994; Hellier 1996; 
de Martino et al. 2008; Brunschweiger et al. 2009). IPs are 
asynchronous systems where the orbital period of the system 
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P or b is larger than the spin period of the white dwarf P sp m- 
The majority of the IPs have P ap i n /Porb < 0.1 , theoreti- 
cally they are predicted to lie in the 0.01-0.6 range (Norton, 
Wynn & Somerscales 2004; Norton et al. 2008; Scaringi et 
al. 2010). IPs are well characterised by high complex ab- 
sorption in their X-ray spectra which extends, where the 
spectra are available in some cases, out to 100 keV. They 
show modulation of the X-ray and/or optical light curves at 
P S pi n and the appearance of the beat period (where 1/Pbeat 
= 1/Pspin - 1/Porb) or its sidebands. Many mCVs are de- 
tected by INTEGRAL IBIS (Barlow et al. 2006; Landi et 
al. 2009), and Swift BAT (Brunschweiger et al. 2009) in the 
20-100 keV ranges and most of them are found to be IPs. 
This is thought to be because the hardest X-ray emission 
comes from the cooling post-shock region of IPs rather than 
Polars, since the strong magnetic field of Polars causes the 
cyclotron cooling to dominate the cooling post-shock region 
suppressing the Bremsstrahlung emission. 

EX Hya was discovered by Kraft (1962) with an orbital 
period of 98 minutes. Later, the 67 minutes rotation period 
of the white dwarf was identified (Vogt et al. 1980), clas- 
sifying the source as an intermediate polar. The WD mass 
is estimated to be in the range of 0.5-0.8 M© as measured 
using optical (Beauerman et al. 2003; Hellier et. al 1987) 
and X-ray (Fujimoto & Ishida 1997) observations. The ele- 
mental abundances are in accordance with solar abundances 
except for iron which is less than solar (0.6 ± 0.2) (Fujimoto 
& Ishida 1997). Using optical observations, the system has 
been found to be at a distance of 64.5 ± 1.2 pc, with a qui- 
escent accretion luminosity of (2.6 ± 0.6) x 10 32 erg s , 
a quiescent accretion rate of (6.2 ± 1.5) x 10 -11 Mq yr~ 
and an inclination of 76.0° - 77.6° (Beuermann et al. 2003). 

EX Hya shows a partial eclipse in its light curves over 
the orbital period. The eclipse was explained by two scenar- 
ios; either occultation of one of the poles for two-pole accre- 
tion or partial occultation of the accretion column by the 
secondary for one-pole accretion (Rosen, Mason & Cordova, 
1988). The source is suggested to have an extended bulge 
on the outer accretion disc and extended/overflowing mate- 
rial originating from the hot spot (Belle et al. 2005). The 
source is observed both in outburst and quiescence. Dur- 
ing the outburst stage, the Doppler tomograms and radial 
velocities indicate that the source shows disc overflow and 
dips in the optical light curve (Mhlahlo et al. 2007a); the 
disc overflow is also present in the quiescent state (Mhlahlo 
et al. 2007b) consistent with the above scenario. The evi- 
dence of a bulge on the accretion disc is also present in the 
X-ray observations. The X-ray binary light curves in the 1- 
20 A range show an extended broad absorption associated 
with the disc which increases as the wavelength increases 
(Hoogerwerf, Brickhouse & Mauche 2005). 

The X-ray spectrum of the post shock emission of EX 
Hya has been previously modelled with multi-temperature 
plasma emission models. One example uses four optically 
thin plasma emission models at four different tempera- 
tures (Allan, Hellier & Beardmore, 1998). Another approach 
uses a cooling- flow model (MKCFLOW) assuming a multi- 
temperature thermal plasma emission with a relatively flat 
emission measure distribution indicating cooling gas from a 
steady-state condition. (Mukai et al. 2003). 

The spin modulation of the X-ray light curve of the 
source can be explained either by the absorption difference 



over the spin phase because of the accretion curtain, or oc- 
cultation of the accretion curtain by the WD. The occulta- 
tion scenario is more favoured because the system has tall 
accretion coloumns (~ 1 Rwd ) (Allan, Hellier & Beard- 
more, 1998). The source also shows spectral variability as 
a function of orbital phase due to absorption (Cropper et. 
al 2002), but the modulations and spectral models have not 
been thoroughly investigated up until this work. 

In this work, the X-ray spectral properties of EX Hya 
are thoroughly investigated. In Section 2, the X-ray observa- 
tions and data preparation are introduced. In Section 3 the 
total spectrum for the 2000 and 2003 data are investigated. 
A composite-model fit for both observations is performed 
and the variation of the spectral parameters between the 
two observations are discussed. In Section 4, the characteris- 
tics of the source spectra in both 2000 and 2003 observations 
are studied over the orbital phase. The trends in the derived 
spectral parameters as a function of the orbital phase as well 
as the differences between the 2000 and 2003 data are com- 
pared. In Section 5, a similar analysis is made, this time as 
a function of the spin phase of EX Hya. In Section 6, the 
results are discussed and finally in Section 7 the summary 
and conclusion of our work is presented. 



2 OBSERVATION AND DATA 

Two different archival observations of EX Hya obtained with 
the XMM-Newton Observatory were used in this study. 
The first one was conducted on 01 July 2000 (OBS ID: 
00111020101) with 30 ksec exposure time and the second 
one on 11 January 2003 (OBS ID:0057740301) with 57 ksec 
exposure time. The XMM-Newton Observatory (Jansen et 
al. 2001) has a pn CCD detector (Striider et al., 2001) and 
two MOS CCD detectors sensitive in the 0.2-15 keV en- 
ergy range at the focus of three European Photon Imaging 
Cameras (EPIC) with a field of view of about half a de- 
gree for each. It also has Reflecting Grating Spectrometers 
(RGS) (den Herder et al. 2001), a high resolution spectrom- 
eter working together with the EPIC detectors sensitive in 
the 0.35-2.5 keV energy range as well as a 30 cm optical mon- 
itor instrument (OM) with an optical/UV telescope (Mason 
et al. 2001). 

We first checked the pile up in the observations using 
epatplotQ tool in sas. The EPIC MOS data were piled up 
in the 2000 observation, hence only EPIC pn data was used 
for both observations in this analysis. The 2000 observation 
was a pointed observation of the source in the small window 
mode for the EPIC pn with a net count rate of 40.1 ± 0.04 
count s _1 in the 0.2-10 keV energy interval. In the 2003 
observation the EPIC pn was used in the extended full frame 
mode and the source is off axis by 14 arc-minutes, the net 
source count rate is 9.7 ± 0.01 count s _1 in the 0.2-10 keV 
energy range. The significantly lower count rate in the 2003 
observation is due to vignetting, caused by the source being 
off-axis. We utilized the CALVIEW tool in SAS and determined 
the percentage of loss in the count rate as 65% which yields 
a count rate of about 13-14 count s _1 based on 40.1 count 
s _1 from 2000. This is in accordance with our count rate 
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derived from the 2003 data. However, we note a difference 
of about 30-40% between the rate calculated from the data 
and our inspection from CALVIEW using the rate of the 2000 
data. This may be a result of extra vignetting (corrected in 
the analysis) or about 30-40% change in the source count 
rate, which is not conclusive at this stage. 

For both observations, the standard pipeline processed 
data was used to perform the spectral and temporal analysis 
using the XMM-Newton Science Analysis Software (sas) 
version 8.0.0. The source light curve and spectrum extrac- 
tion are done using the xmmselect tool within SAS. In the 
2000 observation, a circular photon extraction region was 
centered on the source with a radius 38.8 arcsec, and a same 
size region was used for the background extraction centered 
elsewhere on the image. For the 2003 data, since the source 
was off-axis, an elliptical region with a semi-major axis of 60 
arcsec and a semi-minor axis of 23.3 arcsec, orientated at an 
angle of 320 degrees was used in order to extract the source 
and background events. Data with single- and double-pixel 
events, i.e., patterns 0-4 with Flag=0 options were used. 
In order to perform spectroscopy, we used the PHASECALC 
tool in SAS, and created phase columns in the event files us- 
ing the spin and orbital periods; and spectra were extracted 
from the event files for a phase interval of 0.1 with iterative 
calls to the especget tool in SAS yielding phase- resolved 
spectra, xspec 12.2.1 (Arnaud 1996) was used for further 
analysis. 



3 TOTAL SPECTRUM OF THE SOURCE 

Prior to the fitting procedure, spectra were calculated using 
the extraction criteria described in section 2 in the energy 
range of 0.2-10.0 keV. The resulting spectra were grouped so 
that the minimum counts for each energy bin were 300 for 
the year 2000 observation and 200 for the 2003 observation 
using XSPEC. The fits were conducted in the 0.2-10.0 keV 
energy band. 

The first attempts at fitting the total spectra using 
XSPEC were done by using a plasma emission model that 
is in collisional equilibrium (MEKAL model in xspec) and 
creating a composite model of several MEKAL models to fit 
the spectra at several different temperatures as suggested by 
Allan et al. (1998). Although we included up to four different 
MEKAL model components, the trials yielded fits with re- 
duced x 2 more than 2. Then a different approach was taken 
by using a cooing-flow model (the VMCFLOW model in 
xspec) suggested by Mukai et al. (2003) which again failed 
to yield good fits. Finally, we combined the two methods, 
using a cooling-flow model (VMCFLOW) together with two 
different temperature MEKAL models at around 0.64 keV 
and 1.6 keV. An absorber model for the interstellar extinc- 
tion and a partially covering absorber model (TBABS and 
PCFABS models in xspec) was assumed to account for the 
intervening absorption for both of the observations. 

The fits with the multi-component model yielded fluc- 
tuating residuals around particularly the iron line complex 
at about 6.7 keV. We believe this problem arises from CTI 
effects around particularly strong lines in the EPIC pn spec- 
trum. Thus, we decided to model the lines separately and 
decreased the Fe abundance in the VMCFLOW model to 
0.2 times the solar abundance and added Gaussian emis- 



sion lines to the model at 6.7 keV and 6.9 keV to account 
for the Fea emission line complex. This procedure enabled 
us to individually model the Fe emission lines. We investi- 
gated any effects of the reduction of the Fe abundance in 
modeling the spectra and we detected no other effects. This 
is expected since the VMCFLOW models only the harder 
part of the spectra. Other Gaussian emission lines were in- 
cluded in the fitted composite model at 7.8 keV, 0.58 keV 
and 0.78 keV for the 2000 data and 7.8 keV, 6.4 keV and 
0.58 keV for the 2003 data in order to reduce the high sigma 
deviations in the residuals. The 7.8 keV, 0.58 keV and 0.78 
keV lines correspond to Fe XXV K/3, O VII and O VIII 
transitions respectively according to the CHIANTI atomic 
database. Modeling these lines by increasing the related el- 
ement abundances is not possible, hence Gaussian profiles 
are used instead. The fitted spectra are presented in Figure 
1 and the spectral parameters are displayed in Table 1. 

We calculate unabsorbed fluxes with similar values of 
1.43 x 10~ 10 and 1.24 x 10~ 10 erg cm" 2 s" 1 for 2000 
and 2003 observations, respectively. These flux values yield 
luminosities of 6.7 x 10 31 erg s _1 and 5.8 x 10 31 erg 
s _1 assuming a 64.5 pc source distance. Using the relation 
L=GMM/2R and taking the mass of the white dwarf as 0.5 
Mq and radius as 0.95 Hq (Beuermann et al. 2003), we find 
mass accretion rates of 5.75 x 10 -11 Mq yr _1 and 5.0 x 
10 -11 M Q yr" 1 for 2000 and 2003 data respectively. 

Both spectra show an unusual feature; an emission line 
at around 7.8 keV. To investigate this feature more closely, 
the spectra of the 2000 and 2003 data, were fitted over a 
range of 5-10 keV, with a power law model for the contin- 
uum and three Gaussians to model the emission lines at 
around 6.7 keV, 6.9 keV and 7.8 keV. The resulting param- 
eters of the continuum and the spectral lines are displayed 
in Table 2. We find an emission line centered around 7.8lg' 3 
keV with a line width a = 0.06± keV for 2000, and cen- 
tered around 7.9±° \l with a line width a = 0.5+^;^ for 2003 
observations. We believe this is the Fe XXV K/3 emission 
line. The line widths for the 6.7 keV emission line is similar 
for both observations. In the 2003 data, the line widths for 
6.9 and 7.8 keV lines are significantly higher. 



4 ORBITAL PHASE-RESOLVED 
SPECTROSCOPY 

In this work, a detailed orbital phase-resolved spectroscopy 
of EX Hya was carried out for the first time. For both obser- 
vations, 2000 and 2003, phase columns were created using 
the ephemeris T = 2437699.94179+0. 068233846(4)E given 
by Hcllier & Sproats (1992). Then for each 0.1 phase interval 
from phase to 1, the spectra were extracted as described 
in section 2. The calculated spectra were grouped such that 
each energy bin contained minimum of 100 counts to im- 
prove the statistical quality. The channels below 0.3 and 
above 10.0 keV were ignored during the fitting procedure to 
improve the statistical quality of the spectra. 

To start with, we calculated the ratio of the maximum 
to minimum spectrum of the orbital phase-resolved spectra 
in order to study the effect of absorption/occultation in the 
system in the two different years. The ratio of the normalised 
count rate per energy for phases around the orbital maxi- 
mum (phases between 0.6-0.8) to phases around the orbital 



4 Pekon Y. and Balman, §. 



minimum (phases between 0.9-1.1) are plotted against the 
energy in Figure 2. Both 2000 and 2003 spectral ratios show 
clear spectral variation and a power-law decline as the en- 
ergy increases with ratio values greater than 1 particularly 
at lower energies. There is a more pronounced steep decline 
with increasing energy in the year 2003 showing the exis- 
tence of large scale absorption difference within the system 
from orbital maximum to minimum. 

As in the case of the fits to the total spectra, mod- 
els such as single MEKAL (in xspec), composite multi- 
ple MEKAL or a single VMCFLOW models alone yielded 
large reduced \ 2 values above 2 during the fitting proce- 
dure within XSPEC. Thus, we attempted to fit a combination 
of a MEKAL model with a VMCFLOW using solar abun- 
dances along with a photoelectric absorption model (WABS 
in xspec). However, this approach did not yield good fits 
either. The Fe emission line complex around 6.7 keV still 
yielded fluctuating residuals. Thus, we applied the same pro- 
cedure to reduce the Fe abundance and adding individual 
Gaussian lines to model the 6.7 keV and 6.9 keV iron lines 
when necessary. 

In the analysis of the 2000 data, we used a composite 
model to fit the spectra at different phases comprising of two 
collisional equilibrium plasma emission models, MEKAL1 
and MEKAL2, one cooling-flow plasma emission model, 
VMCFLOW, a photoelectric absorption model, WABS and 
three Gaussians for particular Fe lines. The spectral parame- 
ters from the fit results are displayed in Table 3 and selected 
parameters are plotted against the orbital phase in Figure 3. 
The fitted spectra at orbital minimum (phase 0.9) and max- 
imum (phase 0.3) are shown in Figure 4. For the absorption 
component, a simple cold absorber was enough to fit the 
data, hence a partial covering absorber was not used. The 
column density of neutral Hydrogen (absorption parameter 
(Nh)) shows distinct variation with the orbital phase and 
indicates an anti-correlation with the orbital flux variations 
reaching a maximum at the orbital minimum. This clearly 
shows that orbital variations at X-ray wavelengths are pro- 
duced by absorption on the orbital plane. The LowT (low 
temperature) parameter of the VMCFLOW model shows 
the same anti-correlation with flux and reaches a minimum 
at orbital maximum whereas the HighT (high temperature) 
parameter of the same model shows almost a direct correla- 
tion with the flux which may imply scattering being at work. 
The plasma emission components modeled by the MEKAL 1 
and MEKAL2 models show almost a constant temperature 
of about 0.6-0.8 keV and 1.3-1.7 keV and the normalisation 
varies slightly in phase with the orbital motion. The Gaus- 
sian line at 6.7 keV does not vary with orbital phase apart 
from a slight drop between phases 0.6-0.7 which indicates 
that the region where the line is produced is visible at all 
times during the binary motion. 

In the analysis of the 2003 data, we used a compos- 
ite model comprising two collisional equilibrium plasma 
emission model, MEKAL1 and MEKAL2, one cooling-flow 
plasma emission model, VMCFLOW, a photoelectric ab- 
sorption model, WABS and three Gaussians for particular 
Fe lines. The spectral parameters from the fit results are 
displayed in Table 4 and selected parameters are plotted 
against the orbital phase in Figure 5. The fitted spectra at 
orbital minimum (phase 0) and maximum (phase 0.4) are 
shown in Figure 6. The resulting fits have reduced \ 2 we U 



below 2 for all the phase spectra. We find that some param- 
eters show distinct variations with the orbital phase. The 
Nn parameter for the intrinsic absorption and LowT pa- 
rameter of the VMCFLOW show anti-correlation with the 
flux i.e. they yield higher values at orbital minimum phases 
(phases between 0.9 and 1.3), and lower values at the max- 
imum phases (phases between 0.4 and 0.8). This indicates 
that source emission is being absorbed at particular phases 
(0.9 and 1.3) where there is absorption on the orbital plane 
and the spectrum gets harder in general (softer X-rays are 
absorbed). The HighT parameter shows no distinct varia- 
tion over the orbital phase in the 2003 data. This is obvious 
since harder X-rays will not be affected/absorbed as much 
by the neutral absorption. 

The plasma emission components modeled by the 
MEKAL model in the 2003 data show almost a constant 
temperature of about 0.6-0.8 keV and 1.3-1.7 keV; the nor- 
malisations vary slightly in phase with the flux. These com- 
ponents may be affected by only occultation , i.e., the eclipse 
itself. The width a and normalisation of the 6.7 keV emis- 
sion line shows direct correlation, with low values at orbital 
minimum and high values at orbital maximum. The 6.7 keV 
emission line (Fe XXV line) is present at all phases. For 
the orbital minimum (phases between 0.9 and 1.3), emission 
lines at 6.9 keV (Fe XXVI line) and at 6.4 keV (fluorescent 
Fe line) were necessary for successful fits to the data (See 
Figure 6). Figure 7 shows the necessity to include the 6.4 
and 6.9 keV emission lines in the 2003 orbital spectra at 
minimum phases by fitting 4 MEKAL models with variable 
abundances and focusing on the 6-8 keV regime. The excess 
around 6.4 keV is clear in Figure 7. For the orbital maximum 
(phases between 0.4 and 0.8), the 6.9 and 6.4 keV emission 
lines could not be fitted with free parameters, however, up- 
per limits for normalisations were obtained by fixing all the 
other model parameters. 

The 6.7 keV line widths show no distinct variations in 
the 2000 data, with values around 0.1 keV. On the other 
hand, in the 2003 data, the widths vary in phase with the 
flux; the line is broader at the phases of orbital maximum 
with widths around 0.3 keV, and at phases of orbital mini- 
mum it has similar values to the 2000 data, around 0.1 keV. 
The width of the 6.4 keV line in the 2003 data varies between 
and 0.16 keV. 



5 SPIN PHASE-RESOLVED SPECTROSCOPY 

We have also performed spin phase-resolved spectroscopy in 
a similar fashion to the orbital phase- resolved spectroscopy 
for the 2000 and 2003 data. The ephemeris for the spin pe- 
riod T = 2437699.8914(5) + 0.046546504(9)E given by Hel- 
lier & Sproats (1992) is used for both observations to create 
spin phase-resolved spectra using the method described in 
section 2. xspec was used for the remaining spectral anal- 
ysis, mainly the fitting procedure. The spin phase- resolved 
spectra were grouped so that each energy bin contained a 
minimum of 100 counts and data points below 0.3 keV and 
above 10.0 keV were omitted due to bad statistics. 

The ratio of the spin maximum spectra (phases between 
0.75 and 0.95) to spin minimum spectra (phases between 
0.25 and 0.45) was also studied for the spin phase-resolved 
spectra. For both the 2000 and 2003 data, the maximum 
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to minimum ratios show similar behaviour, thus there is no 
difference in the spin phase-resolved spectra, as far as the 
absorption characteristics are concerned between the two 
different observations of EX Hya. At low energies, the ratios 
have values of about 2 that gradually decrease with increas- 
ing energy down to 1 (see Figure 8). 

Similar to the previous fits to spectra, two MEKAL 
models, a VMCFLOW model and various Gaussians to 
model emission lines were used to fit all the spectra at all 
spin phases (with 0.1 phase increments). Differently from 
the orbital phase-resolved spectroscopy, inclusion of a partial 
covering absorber model (PCFABS in xspec) was necessary 
to achieve good fits in the low energies. 

The 2003 spin phase-resolved spectra were fitted with 
XSPEC models, two MEKAL (collisional equilibrium plasma 
emission model with solar abundances) at different temper- 
atures, VMCFLOW (cooling-flow model), PCFABS and a 
Gaussian for the emission line at 6.7 keV. No other emis- 
sion lines were needed to be included in the fit. Figure 9 
shows the fitted spectra at the spin phase maximum (phase 
0.7) and minimum (phases 0.3). Table 5 displays spectral 
parameters and Figure 10 is the plot of selected spectral pa- 
rameters versus the spin phase for the 2003 data. The lower 
and upper temperature of the cooling flow (i.e. LowT and 
HighT parameters of the VMCFLOW model) show anti- 
correlation with the light curve folded over the spin period. 
The Nh and the covering fraction of the partial covering 
absorber parameter does not show a direct anti-correlation 
with the flux as would be expected from neutral absorption 
producing the X-ray modulations. The spectral parameters 
of the 6.7 keV line do not show any variation over the spin 
phase. As with the orbital phase-resolved spectroscopy, the 
need for emission lines were investigated by fitting the data 
with 4 MEKAL models with varying abundances and close 
inspection over the 6-8 keV range. The excess at 6.4 keV 
does not show any significant pattern as in the orbital phase 
spectra. The 6.4 keV and 6.9 keV emission lines could not 
be fitted with free parameters. However the upper limits for 
the 6.4 and 6.9 keV line normalisations could be calculated, 
with values 1.0 x 10~ 8 and 5.2 x 10~ 7 for 6.9 keV and 2.7 
x 10~ 6 1.4 x 10" 6 for 6.4 keV at minimum flux (phases 0.2 
and 0.3 ); 9.4 x 10~ 10 and 8.5 x 10~ 6 for 6.9 keV and 2.3 x 
10 -6 and 6.0 x 10~ 7 for 6.4 keV at maximum flux (phases 
0.7 and 0.8). This may be consistent with our results since 
we find components and emission lines (like the 6.7 keV line) 
that show widths and normalisations independent of phase 
generally except for orbital phase in 2003. 

The fits to the spin phase-resolved spectra of the 2000 
data show that the 6.7 keV and 6.9 keV emission lines were 
persistent at all phases. In 2000 data, the parameters Nn, 
covering fraction, LowT and HighT have higher values dur- 
ing orbital minimum (phases between 0.1 and 0.4) and lower 
values otherwise showing anti-correlation with the flux. This 
is the generally expected case if absorption is at work creat- 
ing the modulations and the spectrum gets harder towards 
the minimum as the softer photons are absorbed. The spec- 
tral parameters of the 6.7 keV line do not show any variation 
over the spin phase in the 2000 data. The 6.9 keV line is also 
persistent over the entire range of spin phase except that 
the normalisation minimizes around phases 0.4-0.5 which is 
about the spin minimum. Figure 11 shows the fitted spec- 
tra at the spin phase maximum (phase 0.9) and minimum 



(phase 0.5). Table 6 presents the spectral parameters and 
Figure 12 displays the plot of selected spectral parameters 
versus the spin phase. 

The MEKAL emission components do not show any 
significant variation of the temperature parameter over all 
phases, but a slight variation of the normalisation in phase 
with the spin modulation for both 2000 and 2003 observa- 
tions, is seen. 

The widths and normalisations of the emission lines 
show no significant variation with spin phase for both 2000 
and 2003 data, while the width of the 6.7 keV line in 2003 
data is slightly larger than that of 2000. 



6 DISCUSSION 

We have presented the XMM-Newton EPIC pn data of 
the intermediate polar EX Hya. The spectrum of EX Hya is 
quite complex as it needs multi-component emission mod- 
els to fit the data. For the total spectra of both 2000 and 
2003 observations, the composite model is composed of 2 
different MEKAL models with plasma temperatures around 
0.6-0.8 keV and 1.3-1.7 keV, a cooling-flow emission model, 
VMCFLOW with a distribution of 3-33 keV for the 2000 
observation and 8-61 keV for the 2003 observation, and var- 
ious Gaussian emission lines at 6.4 keV, 6.7 keV, 6.9 keV, 
0.58 keV, 0.78 keV and 7.8 keV. The 7.8 keV emission line 
is an unusual feature encountered in the X-ray spectrum 
of CVs. After a closer investigation around the 5-10 keV 
spectra of both observations, the line appeared prominently 
with a centre around 7.8lg'g, line width of 0.2-0.4 keV and 
with a normalisation of 0.0001 photons cm~ 2 s _1 . This line 
may be the Fe XXV K/3 emission line according to the CHI- 
ANTI atomic database. The emission lines at 0.58 and 0.78 
keV correspond to O VII and O VIII transitions respectively. 
These lines appear only in the phase averaged spectra, hence 
may be calibration artifacts. 

The two different collisional equilibrium plasma emis- 
sion components around 0.6-0.8 and 1.3-1.7 keV (modeled 
with MEKAL) are present at all phases (orbital and spin) 
and in the total spectra. This shows that there is an X-ray 
emitting region unaffected by orbital and spin phase which 
may be attributed to an X-ray emission region in the line of 
sight at all phases. We find that the normalisations of the 
MEKAL emission components decrease during the eclipse, 
which indicates that these components originate most likely 
from the post-shock region in the accretion column which 
is too complex to be represented with just a cooling-flow 
plasma emission model alone. 

The folded light curve of the source over the orbital 
period shows two distinct features for both 2000 and 2003 
data. The first feature is the eclipse expected to be seen at 
phases around 0. For the 2000 light curve the eclipse falls at 
phases 0.93 ± 0.05. Fitting a simple sinusoid to the folded 
light curve yields a semi-amplitude variation of 6.4 ± 1.0 
counts s _1 and 14% modulation. For 2003 data the eclipse 
is at phase 1.02 ± 0.05 with semi-amplitude variation of 1.6 
± 0.3 counts s _1 and 20% modulation. The accumulated er- 
rors in orbital phase are 0.000083 for 2000 and 0.000088 for 
2003 data. With these values, it is clear that the two eclipses 
do not overlap, however early works (e.g. Hoogerwerf, Brick- 
house & Mauche 2005) show that the centroid of the eclipse 
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at high energies (i.e. for 10-15 A it is at phase 0.992) falls 
behind the eclipse at low energies (i.e. for 15-20 A it is at 
phase 1.0). Furthermore EUV eclipse was found to be at 
phase 0.97 by Belle et al. (2002) and it was suggested that 
the optical eclipse at phase might be due to the eclipse of 
the hot spot rather than the white dwarf. Hence, the shifting 
of the eclipse may be to the change in geometry of the eclips- 
ing region between two observations. During the eclipse, the 
count rate does not drop to zero, which is consistent with a 
double pole accretion model where the lower accretion col- 
umn is eclipsed but the upper column is not (Belle et al. 
2005) . This is consistent with our results since we find com- 
ponents and emission lines (like the 6.7 keV line) that show 
widths and normalisations independent of any type of phase. 
However, it is inconclusive since the spin-folded light curves 
show a single peak profile at all times. 

The second prominent feature of the folded mean light 
curve over the orbital phase is the broad shallow dip ap- 
pearing in both 2000 and 2003 data. In the 2000 data, the 
dip is at around phases 0.45-0.65 and roughly 1/3 of the 
depth of the eclipse. In the 2003 data, the absorption fea- 
ture is shallower, (about 1/5 of the depth of the eclipse) 
and the location of the dip is around phases 0.35-0.55. This 
broad modulation was well observed over the optical, UV 
and X-ray regime and attributed to the bulge and/or hot 
spot on the disc (e.g. Belle et al. 2005, Mhlahlo et al. 2007, 
Hoogerwerf et al. 2005). The dip gets stronger and moves 
to phases closer to the eclipse as the wavelength increases 
due to the changes in the ionization state of the bulge. In 
the RXTE light curve, the dip is around phase 0.65 (Belle 
et al. 2005) , and in Chandra light curve it is around phase 
0.78 (Hoogerwerf et al. 2005). Both of these observations are 
from May 2000, and the broad dip feature in the July 2000 
XMM-Newton light curve is in accordance with the above 
examples. However the feature in 2003 is not consistent. The 
dip is shallower and at phases further away from the eclipse. 
This may imply that the source has undergone structural 
changes between 2000 and 2003 and the bulge has shifted 
position or grown in size so that modulation is observed at 
earlier phases. Moreover, the phase-resolved spectra show no 
significant change in the absorption for either 2000 or 2003 
observations, hence the modulation in the orbital light curve 
cannot be attributed to neutral absorption on the disc. 

In both 2000 and 2003 data, the neutral absorption and 
low temperature values of the cooling flow increase signifi- 
cantly during the eclipse. This is an indication that cold 
absorption plays an important role in the orbital variation 
of the source. During the eclipse, the absorption increases, 
reducing the emission from the lower end of the flow and 
hence the lower temperature values increase. The change 
in absorption over the orbital motion has been encountered 
before in another CV; a classical nova RR Pic (Pekon & 
Balman 2008) where an absorbing region on the disk in the 
line of sight affects the spectrum. This makes the spectral 
behaviour of CVs over the orbital period an area worth more 
investigation. 

The light curve folded over the spin period of the white 
dwarf shows clear sinusoidal modulation for both 2000 and 
2003 data. A simple sinusoidal fit shows that the spin mini- 
mum is at phase 0.26 ± 0.03, the semi-amplitude variation 
is 2.1 ± 0.1 counts s _1 with 25% modulation for 2003 data. 
For the 2000 data the minimum is at phase 0.36 ± 0.06, the 



semi-amplitude variation is 14.5 ±1.2 counts s 1 with 30% 
modulation. The accumulated errors in phase are 0.00027 for 
2000 and 0.00029 for 2003 data, thus they are much smaller 
than our phase errors. There is a difference in the phase of 
the spin minimum in different years, most likely due to the 
changes of the accretion column/curtain structure. 

The fits show a net increase of absorption coloumn and 
covering fraction at the spin minima in the 2000 data indi- 
cating the role of absorption from the curtains. Moreover, 
the low temperature limits of the plasma increase at the 
phases of minimum flux, making the spectrum harder; since 
absorption reduces the amount of low temperature X-rays. 
In the 2003 data the absorption coloumn and covering frac- 
tion do not show a clear correlation with the flux, but the low 
temperature limits of the plasma modulate inversely with 
the flux. These findings imply that the role of absorption 
from the accretion curtain is dominant for the source. Al- 
though our results do not rule out occultation by the white 
dwarf as proposed by Allan et al. (1998), the findings favour 
the absorption effects on the spin modulation to be more 
dominant. 

X-ray spectra of IPs commonly show the fluorescent 6.4 
keV emission line feature originating from the reprocessing 
of the hard X-rays from the white dwarf surface (e.g. Hellier 
& Mukai 2004, Landi et al. 2008). The 6.4 keV line was en- 
countered in the Chandra EEC spectrum of EX Hya (Hellier 
& Mukai 2004), where the line was weaker than the other 
lines in the Fe K a complex (i.e. 6.7 and 6.9 keV emission 
lines). The total spectrum of the 2003 XMM-Newton ob- 
servations agrees with this, since the normalisation of the 
Gaussian fitted to the 6.4 keV line is less than those of the 
6.7 and 6.9 keV emission lines. The behaviour of this flu- 
orescent line over the orbital cycle is however unusual. In 
the phase-resolved spectra, the line is detected for the or- 
bital and spin phase-resolved spectra of the 2000 observation 
and the spin phase-resolved spectra of the 2003 data. In the 
analysis of the 2003 data over the orbital phases, the line is 
undetected at phases 0.4 to 0.8 (i.e., phases of orbital max- 
ima) but present otherwise (i.e., phases of orbital minima). 
This means that there is a region on the orbital plane re- 
flecting the X-ray emission seen around the eclipse phases 
which is a region on the disc rather than the surface of the 
white dwarf. This reflection phenomenon is detected from a 
CV and an IP for the first time. We strongly suggest that 
the bulge at the accretion impact zone will appear at similar 
binary phases to the eclipse and causes this reflection to oc- 
cur. Such features are detected in dipping Low-mass X-ray 
binaries (Diaz-Trigo et al. 2006), although with significantly 
higher line widths (e.g. a ~ 0.85 keV) than seen in this work 
(<r ~ 0.1 keV). 

It has been shown that a growing number of IPs have a 
soft blackbody emission component with temperatures less 
than 100 eV, (e.g. Evans & Hellier 2007, Anzolin et al. 2008, 
Anzolin et al. 2009, de Martino et al. 2008). EX Hya is 
listed as a "soft-IP" by Evans & Hellier (2007) however not 
with much emphasis, since the spectral fits with a blackbody 
component yield poor results. Our fits with a blackbody 
component also give reduced \ 2 values much greater than 
2, thus we are also unable to list the source as a "soft-IP". 
There is no need for a soft blackbody component to model 
the XMM-Newton data of EX Hya for both observations. 

The cooling flow-plasma emission component yield tern- 
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perature ranges 3-33 keV (for 2000 data) and 8-61 keV (for 
2003 data) in accordance with the high plasma tempera- 
tures obtained from other IPs as expected (e.g. Anzolin et 
al. 2008, Anzolin et al. 2009, de Martino et al. 2008). We 
note that given these observations (e.g., Brunschweiger et 
al. 2009) EX Hya is one of the IPs that show the hardest X- 
ray emission. In previous works with different observations 
the spectra showed different temperature distributions. Such 
measurements include 19.4 keV from SWIFT data (Brun- 
schweiger et al. 2009), 0.54-15.4 keV from ASCA data (Allan 
et al. 1998), < 20 keV from Chandra data (Mukai et al. 2003) 
and were obtained by fitting plasma models to the spectra, 
extending up to ~80 keV ( Chandra HETG data) using mea- 
surements of the broad components of emission lines form 
the pre-shock flow (Luna et al. 2010). 

The count rate ratio plots of the spin maximum to min- 
imum and orbital maximum to minimum gives us an idea 
about the absorber characteristics. The 2000 and 2003 spin 
ratios are similar, starting at a value about 2 and grad- 
ually falling down to 1 with increasing energy. This may 
imply that the absorption due to the accretion curtain has 
similar structure in both 2000 and 2003. However, we find 
that the neutral hydrogen column density is about a fac- 
tor of four larger in the 2000 data compared with that in 
2003. So the absorption in the accretion column/curtain is 
removed to some extent in the 2003 data. The orbital spec- 
tral ratios show a significant diffence in the years 2000 and 
2003. They both show a power-law decay over the XMM- 
Newton energy range, where the 2000 data is flatter and 
the 2003 data shows a prominent decay in the ratio over the 
XMM-Newton energy range. 

So there is a clear difference between 2000 and 2003 
states of the source as far as absorption over the orbital 
cycle is concerned. This most probably due to the changing 
geometry and ionization structure of the absorbing material 
on the binary plane/accretion disc. 

The estimates from the model integrated fluxes yielded 
luminosities of 6.7 x 10 31 and 5.8 x 10 31 erg s _1 and accre- 
tion rates of 5.75 x 10 -11 and 5.0 x 10 -11 M yr" 1 for 2000 
and 2003 data respectively. The X-ray luminosities are lower 
but the accretion rates are in accordance with the previous 
findings of Beuermann et al. (2003) from optical data. We 
find that this region is located around the accretion impact 
zone creating a sizeable absorbing bulge in a similar fashion 
to Low-mass X-ray binary systems. This region has a smaller 
neutral hydrogen column density by at least a factor of four 
in the year 2000 and affects a smaller range of phases, but 
it still exists. 



models around 0.6-0.8 keV and 1.3-1.7 keV and a cooling- 
flow plasma emission model with a distribution of 3 to 33 
keV for the 2000 state and 8 to 61 keV for the 2003 state. 
We find very little change in the spectral parameters of the 
two collisional equilibrium plasma emission models over the 
orbital and spin cycles. In general, the two observations of 
the source reveal spectral variations. For the 2000 data, we 
find neutral hydrogen column densities in a range 0.01-0.2 x 
10 22 cm~ 2 from spin phase-resolved spectroscopy (with cov- 
ering fractions ranging from 0.2-0.5), and 0.01-0.08 x 10 22 
cm~ 2 from orbital phase-resolved spectroscopy (where cov- 
ering fraction is 1.0). For the 2003 data, the neutral hy- 
drogen column densities are in a range 0.03-0. lx 10 22 cm~ 2 
(with covering fractions of 0.6-1.0) in the spin resolved spec- 
troscopy, and 0.02-0.15 x 10 22 cm~ 2 in the orbital phase- 
resolved spectroscopy (with covering fraction of 1.0). The 
highest Ah values always occur either at the orbital or spin 
phase minima. There is a definite difference in the absorber 
geometry for the spin phase-resolved spectroscopy between 
the two epochs. In the year 2000 it shows a covering frac- 
tion of 0.2-0.5 whereas in the year 2003 it shows a covering 
fraction of 0.6-1.0 (both years show that covering fraction 
reduces towards spin maxima). This indicates a geometri- 
cal change in the absorber from 2000 to 2003. The angle 
subtended by the absorber is considerably less in 2000 and 
hence the accretion curtain is reduced in size. Since we at- 
tained a model of only neutral absorption, we can not com- 
ment much on the ionization state changes of the absorber. 
We note that we do not rule out any occultation effect or 
change in the occultation effects. The absorption in the or- 
bital plane changes by a factor of 2 from 2000 to 2003 and the 
absorption effects are more spread over the orbital cycle in 
comparison with the 2000 data. Thus, the bulge at the accre- 
tion impact zone is larger (and optically thicker) in 2003. We 
also note that both the spin and orbital phase-resolved spec- 
troscopy of the 2003 data show a larger variation (3-25 keV) 
for the LowT parameter of the cooling-flow plasma emission 
model as compared to the 2000 state where this parameter 
stays around 6-12 keV over most of the phase- resolved spec- 
troscopy. Thus, the structure of the cooling flow and the 
absorption effects are different in the two different epochs of 
observations. This is not related to any mass accretion rate 
difference, but only to a difference in the accretion column 
temperature distribution structure and absorption. 



7 SUMMARY AND CONCLUSIONS 

The XMM-Newton EPIC pn (0.3-10. keV) spectrum of the 
IP EX Hya together with spin and orbital phase-resolved 
spectroscopy were thoroughly investigated for the first time 
in this work. We have used two different observations of EX 
Hya to look for spectral variations over time. In both obser- 
vations the mass accretion rate and luminosities are compat- 
ible with those of the previous quiescent observations, hence 
the source is in a quiescent state both in 2000 and 2003. 
The X-ray emission is modelled as arising from three differ- 
ent components: two collisional equilibrium plasma emission 



We detect a 6.4 keV fluorescent Fe line at orbital phases 
0.9-1.3 corresponding to the location of orbital minima in 
2003 and it is absent in the orbital phases 0.4 to 0.8. This 
feature can be attributed to reflection from a region on the 
orbital plane most likely the bulge at the accretion impact 
zone and spread around the disk covering a location includ- 
ing the eclipse phases. So the absorber may be viewed in 
other wavelengths at around these phases. We note that 
a 6.4 keV fluorescent Fe line is non-existent in the orbital 
phase-resolved spectroscopy of the 2000 data. This strongly 
indicates an accretion geometry change and a change in the 
bulge size (impact zone) from year 2000 to 2003; the bulge 
is larger and more spread in 2003. 
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Table 1. Spectral parameters of the year 2000 and 2003 observations in the 0.3-10 keV range. Both spec- 
tra were fitted with a composite model of two collisional equilibrium plasma emission models at different 
temperatures (MEKAL), two Gaussians centered at 6.7 keV, and 6.9 keV, a variable cooling-flow plasma 
emission model (VMCFLOW), a partial covering absorber model (PCFABS) and a simple absorption model 
(TBABS). For the year 2000 spectrum, three more Gaussians were added at 0.58 keV, 0.78 keV and 7.8 keV. 
For the year 2003 spectrum three more Gaussians were added at 6.4 keV, 0.58 keV and 7.8 keV. is the 
absorbing neutral hydrogen column density, CovcrFrac is the covering fraction of the absorber. The given 
errors correspond to 2a confidence level for a single parameter. Gaussian line centers were fixed at their best 
fit values and the iron abundance parameter in VMCFLOW was fixed to 0.2 to improve the quality of the 
fits. 



l\/Tnr]ol 
lVlOOCl 


Component 




9finn Tntnl Qr-ior-triiTYi 

zuuu xuLdi ojjecii mil 


iUUO xoiai opecLiLiui 


tbabs 


N H (x 10 22 


atoms/cm 2 ) 


0.0097iro°oo 7 7 


0092+ - 0025 
u.uuy^_ 003g 


pefabs 


N H (x 10 22 


atoms/cm ) 


0.86+°,;!? 


0-i7±g;S| 




CoverFrac 




22+ ' 01 
""^ -0.01 


0.35«i 


MEKAL 1 


kT 




66+ 01 


63+ 01 




norm (x 10" 


" 3 ) 


4 8+°' 2 


2 9+°' 2 


MEKAL2 


kT 




1 58+ ' 05 


1 65+ ' 07 




norm (x 10" 


" 3 ) 


7 8+°' 7 
'•S-0.7 


4 3+0-5 
^-o.s 


Gaussianl (6.7 keV) 


sigma (keV) 




070+ 006 


076+ 009 
U ' U,D -0.010 




norm (x 10" 


- 4 ) 


9 S+° 1 
2 - 8 -0.1 


2 6+°' 2 


Gaussian2 (6.9 keV) 


sigma (keV) 




042+ 016 
u - u4z -0 .020 


0.064±°;°» 




norm (x 10" 


" 5 ) 


6 3+°' 6 


s-sii.i 


Gaussian3 (6.4 keV) 


sigma (keV) 




N/A 


11+"- 04 




norm (x 10" 


" 5 ) 


N/A 


45 +i.i 
^■°-i.o 


Gaussian4 (7.8 keV) 


sigma (keV) 




06+ 06 


o-5il° ; 14 




norm (x 10" 


- 5 ) 


o 7+1-0 
z -< -0.7 


18.211:? 


Gaussian5 (0.58 keV) 


sigma (keV) 




(frozen) 


(frozen) 




norm (x 10" 


-4) 


o n+0.6 


c 9+0.8 
°- z -0.8 


Gaussian6 (0.78 keV) 


sigma (keV) 




(frozen) 


N/A 




norm (x 10" 


- 4 ) 


5 0+°' 6 


N/A 


VMCFLOW 


LowT 




3 4+ 06 


9 0+ 1 ' 2 




HighT 




29.1+11 


S9.2+ 2 ' 4 




Norm (x 10 


-8) 


°- a -0.7 


9 Q" 1 " -' 2 




xl 




1.70 (956 d.o.f.) 


1.66 (903 d.o.f.) 



Table 2. Parameters for the common Gaussian emission lines in 2000 and 2003 spectra. 
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Table 3. Spectral parameters derived for the year 2000 observation at each orbital phase of 0.1 in the 0.3-10 
keV range. All the spectra were fitted with a composite model of two collisional equilibrium plasma emission 
models at different temperatures (MEKAL), 2 Gaussians centered at 6.7 keV and 6.9 keV, a variable cooling- 
flow plasma emission model (VMCFLOW) and photoelectric absorption of HI (WABS). The given errors 
correspond to 2a confidence level for a single parameter. Gaussian line centers were fixed at their best fit 
values and the iron abundance parameter in VMCFLOW was fixed to 0.2 to improve the quality of the fits. 
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Table 4. Spectral parameters derived for the year 2003 observation at each orbital phase of 0.1 in the 
0.3-10 keV range. All the spectra were fitted with a composite model of two collisional equilibrium plasma 
emission models at different temperatures (MEKAL), 2 Gaussians centered at 6.7 keV and 6.9 keV, a 
variable cooling- flow plasma emission model (VMCFLOW) and photoelectric absorption of HI (WABS). 
In some phases another Gaussian at 6.4 keV was added to the model. The given errors correspond to 2tr 
confidence level for a single parameter. Gaussian line centers were fixed at their best fit values and the iron 
abundance parameter in VMCFLOW was fixed to 0.2 to improve the quality of the fits. 
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Table 5. Spectral parameters derived for the year 2003 observation at each spin phase of 0.1 in the 0.3-10 
keV range. All the spectra were fitted with a composite model of two collisional equilibrium plasma emission 
models at different temperatures (MEKAL), a Gaussian centered at 6.7 keV, a variable cooling-flow plasma 
emission model (VMCFLOW) and a partial covering absorber model (PCFABS). The given errors correspond 
to 2<j confidence level for a single parameter. Gaussian line centers were fixed at their best fit values and the 
iron abundance parameter in VMCFLOW was fixed to 0.2 to improve the quality of the fits. 
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Table 6. Spectral parameters derived for the year 2000 observation at each spin phase of 0.1 in the 
0.3-10 keV range. All the spectra were fitted with a composite model of two collisional equilibrium plasma 
emission models at different temperatures (MEKAL), 2 Gaussians centered at 6.7 keV and 6.9 keV, a variable 
cooling- flow plasma emission model (VMCFLOW) and a partial covering absorber model (PCFABS). The 
given errors correspond to 2cr confidence level for a single parameter. Gaussian line centers were fixed at their 
best fit values and the iron abundance parameter in VMCFLOW was fixed to 0.2 to improve the quality of 
the fits. 
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Figure 1. The top panel shows the total EPIC pn spectrum of the year 2000 data (see Table 1 for details). The bottom panel shows 
the total EPIC pn spectrum of the year 2003 data (see Table 1 for details). 
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Figure 2. The left-hand panel shows the ratio of the count rate spectrum of the orbital maximum to the spectrum of the orbital minimum 
in the year 2000 data. The right-hand panel shows the ratio of the count rate spectrum of the orbital maximum to the spectrum of the 
orbital minimum in the year 2003 data. 
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Figure 3. The light curve of the year 2000 data folded over the orbital period of 98 minutes using the ephemcris T = 
2437699. 94179+0. 068233846(4)E (top panel) and plots of the spectral parameters derived from the orbital phase- resolved spectroscopy 
over the same orbital phase range. The CGauss an d norm Gauss arc parameters for the 6.7 keV line. The lower boundary errors were used 
in plotting the highT parameter. 
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Figure 4. The top panel shows the composite model (WABS*(MEKAL+MEKAL+GAUSS+GAUSS+VMCFLOW)) fitted to the spec- 
trum of orbital maximum (phase 0.3) in the year 2000 data. The bottom panel shows the same composite model fitted to the spectrum 
of the orbital minimum (phase 0.9) in the year 2000 data. The crosses show the data with error bars, solid lines show the composite 
model, the dotted lines show the individual models and the panels underneath show the residuals in standard deviations. 
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Figure 5. The light curve of the year 2003 data folded over the orbital period of 98 minutes using the cphcmcris T = 
2437699. 94179+0. 068233846(4)E (top panel) and plots of the spectral parameters derived from the orbital phase- resolved spectroscopy 
over the same orbital phase range. The CGauss an d norm Gauss arc parameters for the 6.7 keV line. The lower boundary errors were used 
in plotting the highT parameter. 
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Figure 6. The top panel shows the composite model (PCFABS*(MEKAL+MEKAL+GAUSS+VMCFLOW)) fitted to the spectrum of 
orbital maximum (phase 0.4) in the year 2003 data. Note that only 6.7 keV line is present. The bottom panel shows the same composite 
model fitted to the spectrum of the orbital minimum (phase 0) in the year 2003 data. Note that 6.4 keV, 6.7 keV and 6.9 keV Fe emission 
lines are present. The crosses show the data with error bars, solid lines show the composite model, the dotted lines show the individual 
models and the panels underneath show the residuals in standard deviations. 
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Figure 7. The orbital phase- resolved spectra of 2003 data between phases 0.9 to 1.3 from left to right and down, respectively. The 
energy range is between 6 and 8 keV to focus on the emission lines. The spectra were fitted with 4 MEKAL models to illustrate the 
excess at 6.4 keV. 
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Figure 8. The left-hand panel shows the ratio of the count rate spectrum of the spin maximum phase to the spectrum of the spin 
minimum phase in the year 2000 data. The right-hand panel shows the ratio of the count rate spectrum of the spin maximum to the 
spectrum of the spin minimum in the year 2003 data. 
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Figure 9. The top panel displays the composite model (PCFABS*(MEKAL+MEKAL+GAUSS+VMCFLOW)) fitted to the spectrum 
of the spin maximum (phase 0.7) in the year 2003 data. The bottom panel displays the same composite model fitted to the spectrum of 
the spin minimum (phase 0.3) in the year 2003 data. The crosses show the data with error bars, solid lines show the composite model, 
the dotted lines show the individual models and the panels underneath show the residuals in standard deviations. 
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Figure 10. The light curve of the year 2003 data folded over the spin period of 67 minutes using the ephemeris T = 2437699.8914(5) 
+ 0.046546504(9)E (top panel) and plots of the spectral parameters derived from the spin phase-resolved spectroscopy over the same 
spin phase range. The <TGauss an d norm Gauss are parameters for the 6.7 keV line. The lower boundary errors were used in plotting the 
highT parameter. 
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Figure 11. The top panel displays the composite model (PCFABS*(MEKAL+MEKAL+GAUSS+GAUSS+VMCFLOW)) fitted to the 
spectrum of the spin maximum (phase 0.9) in the year 2000 data. The 6.7 keV and 6.9 keV emission lines arc explicitly shown. The 
bottom panel displays the same composite model fitted to the spectrum of the spin minimum (phase 0.5) in the year 2000 data. The 6.7 
keV and 6.9 keV emission lines are explicitly shown. The crosses show the data with error bars, solid lines show the composite model, 
the dotted lines show the individual models and the panels underneath show the residuals in standard deviations. 
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Figure 12. The light curve of the year 2000 data folded over the spin period of 67 minutes using the ephemeris T = 2437699.8914(5) 
+ 0.046546504(9)E (top panel) and plots of the spectral parameters derived from the spin phase-resolved spectroscopy over the same 
spin phase range. The <TGauss an d norm Gauss are parameters for the 6.7 keV line. The lower boundary errors were used in plotting the 
highT parameter. 



